\(\int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx\) [33]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (warning: unable to verify)
   Maple [B] (warning: unable to verify)
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 39, antiderivative size = 252 \[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\frac {g E\left (\arcsin \left (\frac {\cos (e+f x)}{1+\sin (e+f x)}\right )|-\frac {a-b}{a+b}\right ) \sqrt {\frac {\sin (e+f x)}{1+\sin (e+f x)}} \sqrt {a+b \sin (e+f x)}}{(a-b) c f \sqrt {g \sin (e+f x)} \sqrt {\frac {a+b \sin (e+f x)}{(a+b) (1+\sin (e+f x))}}}-\frac {2 \sqrt {a+b} \sqrt {g} \sqrt {\frac {a (1-\csc (e+f x))}{a+b}} \sqrt {\frac {a (1+\csc (e+f x))}{a-b}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {g} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {g \sin (e+f x)}}\right ),-\frac {a+b}{a-b}\right ) \tan (e+f x)}{(a-b) c f} \]

[Out]

g*EllipticE(cos(f*x+e)/(1+sin(f*x+e)),((-a+b)/(a+b))^(1/2))*(sin(f*x+e)/(1+sin(f*x+e)))^(1/2)*(a+b*sin(f*x+e))
^(1/2)/(a-b)/c/f/(g*sin(f*x+e))^(1/2)/((a+b*sin(f*x+e))/(a+b)/(1+sin(f*x+e)))^(1/2)-2*EllipticF(g^(1/2)*(a+b*s
in(f*x+e))^(1/2)/(a+b)^(1/2)/(g*sin(f*x+e))^(1/2),((-a-b)/(a-b))^(1/2))*(a+b)^(1/2)*g^(1/2)*(a*(1-csc(f*x+e))/
(a+b))^(1/2)*(a*(1+csc(f*x+e))/(a-b))^(1/2)*tan(f*x+e)/(a-b)/c/f

Rubi [A] (verified)

Time = 0.35 (sec) , antiderivative size = 252, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.077, Rules used = {3015, 2895, 3011} \[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\frac {g \sqrt {\frac {\sin (e+f x)}{\sin (e+f x)+1}} \sqrt {a+b \sin (e+f x)} E\left (\arcsin \left (\frac {\cos (e+f x)}{\sin (e+f x)+1}\right )|-\frac {a-b}{a+b}\right )}{c f (a-b) \sqrt {g \sin (e+f x)} \sqrt {\frac {a+b \sin (e+f x)}{(a+b) (\sin (e+f x)+1)}}}-\frac {2 \sqrt {g} \sqrt {a+b} \tan (e+f x) \sqrt {\frac {a (1-\csc (e+f x))}{a+b}} \sqrt {\frac {a (\csc (e+f x)+1)}{a-b}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {g} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {g \sin (e+f x)}}\right ),-\frac {a+b}{a-b}\right )}{c f (a-b)} \]

[In]

Int[Sqrt[g*Sin[e + f*x]]/(Sqrt[a + b*Sin[e + f*x]]*(c + c*Sin[e + f*x])),x]

[Out]

(g*EllipticE[ArcSin[Cos[e + f*x]/(1 + Sin[e + f*x])], -((a - b)/(a + b))]*Sqrt[Sin[e + f*x]/(1 + Sin[e + f*x])
]*Sqrt[a + b*Sin[e + f*x]])/((a - b)*c*f*Sqrt[g*Sin[e + f*x]]*Sqrt[(a + b*Sin[e + f*x])/((a + b)*(1 + Sin[e +
f*x]))]) - (2*Sqrt[a + b]*Sqrt[g]*Sqrt[(a*(1 - Csc[e + f*x]))/(a + b)]*Sqrt[(a*(1 + Csc[e + f*x]))/(a - b)]*El
lipticF[ArcSin[(Sqrt[g]*Sqrt[a + b*Sin[e + f*x]])/(Sqrt[a + b]*Sqrt[g*Sin[e + f*x]])], -((a + b)/(a - b))]*Tan
[e + f*x])/((a - b)*c*f)

Rule 2895

Int[1/(Sqrt[(d_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Simp[-2*(
Tan[e + f*x]/(a*f))*Rt[(a + b)/d, 2]*Sqrt[a*((1 - Csc[e + f*x])/(a + b))]*Sqrt[a*((1 + Csc[e + f*x])/(a - b))]
*EllipticF[ArcSin[Sqrt[a + b*Sin[e + f*x]]/Sqrt[d*Sin[e + f*x]]/Rt[(a + b)/d, 2]], -(a + b)/(a - b)], x] /; Fr
eeQ[{a, b, d, e, f}, x] && NeQ[a^2 - b^2, 0] && PosQ[(a + b)/d]

Rule 3011

Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/(Sqrt[(g_.)*sin[(e_.) + (f_.)*(x_)]]*((c_) + (d_.)*sin[(e_.) +
(f_.)*(x_)])), x_Symbol] :> Simp[(-Sqrt[a + b*Sin[e + f*x]])*(Sqrt[d*(Sin[e + f*x]/(c + d*Sin[e + f*x]))]/(d*f
*Sqrt[g*Sin[e + f*x]]*Sqrt[c^2*((a + b*Sin[e + f*x])/((a*c + b*d)*(c + d*Sin[e + f*x])))]))*EllipticE[ArcSin[c
*(Cos[e + f*x]/(c + d*Sin[e + f*x]))], (b*c - a*d)/(b*c + a*d)], x] /; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[
b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && EqQ[c^2 - d^2, 0]

Rule 3015

Int[Sqrt[(g_.)*sin[(e_.) + (f_.)*(x_)]]/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*((c_) + (d_.)*sin[(e_.) +
(f_.)*(x_)])), x_Symbol] :> Dist[(-a)*(g/(b*c - a*d)), Int[1/(Sqrt[g*Sin[e + f*x]]*Sqrt[a + b*Sin[e + f*x]]),
x], x] + Dist[c*(g/(b*c - a*d)), Int[Sqrt[a + b*Sin[e + f*x]]/(Sqrt[g*Sin[e + f*x]]*(c + d*Sin[e + f*x])), x],
 x] /; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b*c - a*d, 0] && (EqQ[a^2 - b^2, 0] || EqQ[c^2 - d^2, 0])

Rubi steps \begin{align*} \text {integral}& = -\frac {g \int \frac {\sqrt {a+b \sin (e+f x)}}{\sqrt {g \sin (e+f x)} (c+c \sin (e+f x))} \, dx}{a-b}+\frac {(a g) \int \frac {1}{\sqrt {g \sin (e+f x)} \sqrt {a+b \sin (e+f x)}} \, dx}{(a-b) c} \\ & = \frac {g E\left (\arcsin \left (\frac {\cos (e+f x)}{1+\sin (e+f x)}\right )|-\frac {a-b}{a+b}\right ) \sqrt {\frac {\sin (e+f x)}{1+\sin (e+f x)}} \sqrt {a+b \sin (e+f x)}}{(a-b) c f \sqrt {g \sin (e+f x)} \sqrt {\frac {a+b \sin (e+f x)}{(a+b) (1+\sin (e+f x))}}}-\frac {2 \sqrt {a+b} \sqrt {g} \sqrt {\frac {a (1-\csc (e+f x))}{a+b}} \sqrt {\frac {a (1+\csc (e+f x))}{a-b}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {g} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {g \sin (e+f x)}}\right ),-\frac {a+b}{a-b}\right ) \tan (e+f x)}{(a-b) c f} \\ \end{align*}

Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(4464\) vs. \(2(252)=504\).

Time = 28.85 (sec) , antiderivative size = 4464, normalized size of antiderivative = 17.71 \[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\text {Result too large to show} \]

[In]

Integrate[Sqrt[g*Sin[e + f*x]]/(Sqrt[a + b*Sin[e + f*x]]*(c + c*Sin[e + f*x])),x]

[Out]

(2*Sin[(e + f*x)/2]*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])*Sqrt[g*Sin[e + f*x]]*Sqrt[a + b*Sin[e + f*x]])/((a -
 b)*f*(c + c*Sin[e + f*x])) + (Cot[(e + f*x)/2]*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^2*Sqrt[g*Sin[e + f*x]]*S
qrt[a + b*Sin[e + f*x]]*(-1/2*(a*Sqrt[Sin[e + f*x]])/((a - b)*Sqrt[a + b*Sin[e + f*x]]) - (b*Sqrt[Sin[e + f*x]
])/(2*(a - b)*Sqrt[a + b*Sin[e + f*x]]) - (b*Cos[(3*(e + f*x))/2]*Sec[(e + f*x)/2]*Sqrt[Sin[e + f*x]])/(2*(a -
 b)*Sqrt[a + b*Sin[e + f*x]]) - (b*Sec[(e + f*x)/2]*Sqrt[Sin[e + f*x]]*Sin[(3*(e + f*x))/2])/(2*(a - b)*Sqrt[a
 + b*Sin[e + f*x]]) + (a*Sqrt[Sin[e + f*x]]*Tan[(e + f*x)/2])/(2*(a - b)*Sqrt[a + b*Sin[e + f*x]]) - (b*Sqrt[S
in[e + f*x]]*Tan[(e + f*x)/2])/(2*(a - b)*Sqrt[a + b*Sin[e + f*x]]))*(-2*Tan[(e + f*x)/2]*(1 + Tan[(e + f*x)/2
]) + (2*Sqrt[-a^2 + b^2]*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b^2)]*(-(EllipticE[ArcSin[Sqr
t[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(-b + Sqrt[-a^
2 + b^2])]*Tan[(e + f*x)/2]) + EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b
^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]
*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2]))]))/((a + b*Sin[e + f*x])*Sqrt[(a*Tan[(e + f*x)/2])/(-b +
Sqrt[-a^2 + b^2])])))/(2*(a - b)*f*(c + c*Sin[e + f*x])*((b*Cos[e + f*x]*Cot[(e + f*x)/2]*Sqrt[Sin[e + f*x]]*(
-2*Tan[(e + f*x)/2]*(1 + Tan[(e + f*x)/2]) + (2*Sqrt[-a^2 + b^2]*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x
]))/(a^2 - b^2)]*(-(EllipticE[ArcSin[Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[
2]], (2*Sqrt[-a^2 + b^2])/(-b + Sqrt[-a^2 + b^2])]*Tan[(e + f*x)/2]) + EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 +
b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sqrt[(a*Ta
n[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2]))]))/((a + b*Sin[e
+ f*x])*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])])))/(4*(a - b)*Sqrt[a + b*Sin[e + f*x]]) + (Cos[e +
f*x]*Cot[(e + f*x)/2]*Sqrt[a + b*Sin[e + f*x]]*(-2*Tan[(e + f*x)/2]*(1 + Tan[(e + f*x)/2]) + (2*Sqrt[-a^2 + b^
2]*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b^2)]*(-(EllipticE[ArcSin[Sqrt[(-b + Sqrt[-a^2 + b^
2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(-b + Sqrt[-a^2 + b^2])]*Tan[(e + f*
x)/2]) + EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt
[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]*Sqrt[-((a*Tan[(e + f*
x)/2])/(b + Sqrt[-a^2 + b^2]))]))/((a + b*Sin[e + f*x])*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])])))/
(4*(a - b)*Sqrt[Sin[e + f*x]]) - (Csc[(e + f*x)/2]^2*Sqrt[Sin[e + f*x]]*Sqrt[a + b*Sin[e + f*x]]*(-2*Tan[(e +
f*x)/2]*(1 + Tan[(e + f*x)/2]) + (2*Sqrt[-a^2 + b^2]*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b
^2)]*(-(EllipticE[ArcSin[Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt
[-a^2 + b^2])/(-b + Sqrt[-a^2 + b^2])]*Tan[(e + f*x)/2]) + EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan
[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sqrt[(a*Tan[(e + f*x)/
2])/(-b + Sqrt[-a^2 + b^2])]*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2]))]))/((a + b*Sin[e + f*x])*Sqrt
[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])])))/(4*(a - b)) + (Cot[(e + f*x)/2]*Sqrt[Sin[e + f*x]]*Sqrt[a +
b*Sin[e + f*x]]*(-(Sec[(e + f*x)/2]^2*Tan[(e + f*x)/2]) - Sec[(e + f*x)/2]^2*(1 + Tan[(e + f*x)/2]) - (a*Sqrt[
-a^2 + b^2]*Sec[(e + f*x)/2]^2*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b^2)]*(-(EllipticE[ArcS
in[Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(-b + Sq
rt[-a^2 + b^2])]*Tan[(e + f*x)/2]) + EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a
^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 +
b^2])]*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2]))]))/(2*(-b + Sqrt[-a^2 + b^2])*(a + b*Sin[e + f*x])*
((a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2]))^(3/2)) - (2*b*Sqrt[-a^2 + b^2]*Cos[e + f*x]*Sqrt[(a*Sec[(e + f*
x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b^2)]*(-(EllipticE[ArcSin[Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2
])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(-b + Sqrt[-a^2 + b^2])]*Tan[(e + f*x)/2]) + EllipticF[Arc
Sin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqr
t[-a^2 + b^2])]*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2
+ b^2]))]))/((a + b*Sin[e + f*x])^2*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]) + (Sqrt[-a^2 + b^2]*((
a*b*Cos[e + f*x]*Sec[(e + f*x)/2]^2)/(a^2 - b^2) + (a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x])*Tan[(e + f*x)/2]
)/(a^2 - b^2))*(-(EllipticE[ArcSin[Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]
], (2*Sqrt[-a^2 + b^2])/(-b + Sqrt[-a^2 + b^2])]*Tan[(e + f*x)/2]) + EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^
2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sqrt[(a*Tan[
(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2]))]))/((a + b*Sin[e +
f*x])*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b^2)]*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2
+ b^2])]) + (2*Sqrt[-a^2 + b^2]*Sqrt[(a*Sec[(e + f*x)/2]^2*(a + b*Sin[e + f*x]))/(a^2 - b^2)]*(-1/2*(EllipticE
[ArcSin[Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(-b
 + Sqrt[-a^2 + b^2])]*Sec[(e + f*x)/2]^2) - (a*EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2
])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*Sqrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sec[(e + f*x)/2]^2*Sqrt[(a*Tan[(e
+ f*x)/2])/(-b + Sqrt[-a^2 + b^2])])/(4*(b + Sqrt[-a^2 + b^2])*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^
2]))]) + (a*EllipticF[ArcSin[Sqrt[(b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]/Sqrt[2]], (2*S
qrt[-a^2 + b^2])/(b + Sqrt[-a^2 + b^2])]*Sec[(e + f*x)/2]^2*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2])
)])/(4*(-b + Sqrt[-a^2 + b^2])*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])]) + (a*Sec[(e + f*x)/2]^2*Tan
[(e + f*x)/2]*Sqrt[1 - (-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])])/(4*Sqrt[2]*Sqrt[
-a^2 + b^2]*Sqrt[(-b + Sqrt[-a^2 + b^2] - a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]*Sqrt[1 - (-b + Sqrt[-a^2 + b^2
] - a*Tan[(e + f*x)/2])/(2*Sqrt[-a^2 + b^2])]) + (a*Sec[(e + f*x)/2]^2*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a
^2 + b^2])]*Sqrt[-((a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2]))])/(4*Sqrt[2]*Sqrt[-a^2 + b^2]*Sqrt[(b + Sqrt[-
a^2 + b^2] + a*Tan[(e + f*x)/2])/Sqrt[-a^2 + b^2]]*Sqrt[1 - (b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/(2*Sqr
t[-a^2 + b^2])]*Sqrt[1 - (b + Sqrt[-a^2 + b^2] + a*Tan[(e + f*x)/2])/(b + Sqrt[-a^2 + b^2])])))/((a + b*Sin[e
+ f*x])*Sqrt[(a*Tan[(e + f*x)/2])/(-b + Sqrt[-a^2 + b^2])])))/(2*(a - b))))

Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(3056\) vs. \(2(232)=464\).

Time = 2.48 (sec) , antiderivative size = 3057, normalized size of antiderivative = 12.13

method result size
default \(\text {Expression too large to display}\) \(3057\)

[In]

int((g*sin(f*x+e))^(1/2)/(c+c*sin(f*x+e))/(a+b*sin(f*x+e))^(1/2),x,method=_RETURNVERBOSE)

[Out]

1/2/c/f*(g/((1-cos(f*x+e))^2*csc(f*x+e)^2+1)*(csc(f*x+e)-cot(f*x+e)))^(1/2)*((1-cos(f*x+e))^2*csc(f*x+e)^2+1)*
((a*(1-cos(f*x+e))^2*csc(f*x+e)^2+2*b*(csc(f*x+e)-cot(f*x+e))+a)/((1-cos(f*x+e))^2*csc(f*x+e)^2+1))^(1/2)*((1/
(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)
-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticF((1/(
b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-
a^2+b^2)^(1/2))^(1/2))*(-a^2+b^2)^(1/2)*2^(1/2)*a*(csc(f*x+e)-cot(f*x+e))-2*(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*
x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b
))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticE((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x
+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*(-a^2+b^
2)^(1/2)*2^(1/2)*b*(csc(f*x+e)-cot(f*x+e))-(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)
+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))
*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticF((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+
b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*2^(1/2)*a^2*(csc(f*x+e)-cot(f*x+e))+(1/(b
+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-c
ot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticF((1/(b+
(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^
2+b^2)^(1/2))^(1/2))*2^(1/2)*a*b*(csc(f*x+e)-cot(f*x+e))+2*(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+
(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(
-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticE((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(
-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*2^(1/2)*a^2*(csc(f*x+e)-c
ot(f*x+e))-2*(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)
*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2
)*EllipticE((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^
2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*2^(1/2)*b^2*(csc(f*x+e)-cot(f*x+e))+(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x
+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b)
)^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticF((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+
e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*(-a^2+b^2
)^(1/2)*2^(1/2)*a-2*(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2
)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)
))^(1/2)*EllipticE((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*(
(b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*(-a^2+b^2)^(1/2)*2^(1/2)*b-(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e
)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*2^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1
/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticF((1/(b+(-a^2+b^2)^(1/2))*(a*(cs
c(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*a^2
+(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*
x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticF(
(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2)
)/(-a^2+b^2)^(1/2))^(1/2))*2^(1/2)*a*b+2*(1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b
))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(
csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticE((1/(b+(-a^2+b^2)^(1/2))*(a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b)
)^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*2^(1/2)*a^2-2*(1/(b+(-a^2+b^2)^(1/2))*(a*(c
sc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2)*(1/(-a^2+b^2)^(1/2)*(-a*(csc(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1
/2)-b))^(1/2)*(-a/(b+(-a^2+b^2)^(1/2))*(csc(f*x+e)-cot(f*x+e)))^(1/2)*EllipticE((1/(b+(-a^2+b^2)^(1/2))*(a*(cs
c(f*x+e)-cot(f*x+e))+(-a^2+b^2)^(1/2)+b))^(1/2),1/2*2^(1/2)*((b+(-a^2+b^2)^(1/2))/(-a^2+b^2)^(1/2))^(1/2))*2^(
1/2)*b^2-2*csc(f*x+e)^3*a^2*(1-cos(f*x+e))^3-4*csc(f*x+e)^2*a*b*(1-cos(f*x+e))^2-2*a^2*(csc(f*x+e)-cot(f*x+e))
)/(-cot(f*x+e)+csc(f*x+e)+1)/(a*(1-cos(f*x+e))^2*csc(f*x+e)^2+2*b*(csc(f*x+e)-cot(f*x+e))+a)/(1-cos(f*x+e))*si
n(f*x+e)*2^(1/2)/(a-b)/a

Fricas [F]

\[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\int { \frac {\sqrt {g \sin \left (f x + e\right )}}{\sqrt {b \sin \left (f x + e\right ) + a} {\left (c \sin \left (f x + e\right ) + c\right )}} \,d x } \]

[In]

integrate((g*sin(f*x+e))^(1/2)/(c+c*sin(f*x+e))/(a+b*sin(f*x+e))^(1/2),x, algorithm="fricas")

[Out]

integral(-sqrt(b*sin(f*x + e) + a)*sqrt(g*sin(f*x + e))/(b*c*cos(f*x + e)^2 - (a + b)*c*sin(f*x + e) - (a + b)
*c), x)

Sympy [F]

\[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\frac {\int \frac {\sqrt {g \sin {\left (e + f x \right )}}}{\sqrt {a + b \sin {\left (e + f x \right )}} \sin {\left (e + f x \right )} + \sqrt {a + b \sin {\left (e + f x \right )}}}\, dx}{c} \]

[In]

integrate((g*sin(f*x+e))**(1/2)/(c+c*sin(f*x+e))/(a+b*sin(f*x+e))**(1/2),x)

[Out]

Integral(sqrt(g*sin(e + f*x))/(sqrt(a + b*sin(e + f*x))*sin(e + f*x) + sqrt(a + b*sin(e + f*x))), x)/c

Maxima [F]

\[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\int { \frac {\sqrt {g \sin \left (f x + e\right )}}{\sqrt {b \sin \left (f x + e\right ) + a} {\left (c \sin \left (f x + e\right ) + c\right )}} \,d x } \]

[In]

integrate((g*sin(f*x+e))^(1/2)/(c+c*sin(f*x+e))/(a+b*sin(f*x+e))^(1/2),x, algorithm="maxima")

[Out]

integrate(sqrt(g*sin(f*x + e))/(sqrt(b*sin(f*x + e) + a)*(c*sin(f*x + e) + c)), x)

Giac [F]

\[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\int { \frac {\sqrt {g \sin \left (f x + e\right )}}{\sqrt {b \sin \left (f x + e\right ) + a} {\left (c \sin \left (f x + e\right ) + c\right )}} \,d x } \]

[In]

integrate((g*sin(f*x+e))^(1/2)/(c+c*sin(f*x+e))/(a+b*sin(f*x+e))^(1/2),x, algorithm="giac")

[Out]

integrate(sqrt(g*sin(f*x + e))/(sqrt(b*sin(f*x + e) + a)*(c*sin(f*x + e) + c)), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {g \sin (e+f x)}}{\sqrt {a+b \sin (e+f x)} (c+c \sin (e+f x))} \, dx=\int \frac {\sqrt {g\,\sin \left (e+f\,x\right )}}{\sqrt {a+b\,\sin \left (e+f\,x\right )}\,\left (c+c\,\sin \left (e+f\,x\right )\right )} \,d x \]

[In]

int((g*sin(e + f*x))^(1/2)/((a + b*sin(e + f*x))^(1/2)*(c + c*sin(e + f*x))),x)

[Out]

int((g*sin(e + f*x))^(1/2)/((a + b*sin(e + f*x))^(1/2)*(c + c*sin(e + f*x))), x)